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ABSTRACT 

We investigate the gravitational lensing properties of lines of sight containing multiple cluster- 
scale hales, motivated by their ability to lens very high-redshift {z ~ 10) sources into detectability. 
We control for the total mass along the line of sight, isolating the effects of distributing the mass 
among multiple halos and of varying the physical properties of the halos. Our results show that 
multiple-halo lines of sight can increase the magnified source-plane region compared to the single 
cluster lenses typically targeted for lensing studies, and thus are generally better fields for detecting 
very high-redshift sources. The configurations that result in optimal lensing cross sections benefit 
from interactions between the lens potentials of the halos when they overlap somewhat on the sky, 
creating regions of high magnification in the source plane not present when the halos are considered 
individually. The effect of these interactions on the lensing cross section can even be comparable to 
changing the total mass of the lens from IQ^^Mq to 3 x IQ^^Mq. The gain in lensing cross section 
increases as the mass is split into more halos, provided that the lens potentials are projected close 
enough to interact with each other. A nonzero projected halo angular separation, equal halo mass 
ratio, and high projected halo concentration are the best mass configurations, whereas projected halo 
ellipticity, halo triaxiality, and the relative orientations of the halos are less important. Such high 
mass, multiple-halo lines of sight exist in the SDSS. 
Subject headings: gravitational lensing: strong — galaxies: high-redshift 



1. INTRODUCTION 

Detecting very high-redshift galaxies and determining 
their abundance and star-formation properties is one of 
the foremost challenges in observational cosmology to- 
day. These early galaxies contain the first generation of 
stars in the universe, and their UV emission is thought to 
contribute to the reionization of neutral hydrogen dur- 
ing this epoch. However, observing these faint objects is 
extremely difficult. Past observational studies of z > 7 
galaxies in blank fields using the Lyman-break dropout 
technique (e.g.. iBouwe ns et al. 2004, 2008, 2010a.'^ 



20041: iSchae rer fc Pelld 2005; Richard et al. 2006, 200 
Stark et al.i.2007: .Bradley et al...2008 2012; .Zheng et al 
20091: ILaporte et al.ll2011l: IHah et al.ll20Ta . While lens 
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CasteUano et al."2010'; 'Ouchi et al."2009'; 'Wilkins et alj 
2010, 2011a.b; Capak et al. 2011; Lorenzoni et al. 20lff 
Trenti et al.ll201lL 120121 : lYan et all 1201 lallbl : iHsieh et ail 
2012f). narrow-b a nd searches for Lyman-g emission (e.g. , 
Ota et all [20081: ISobral eFall [200l IHibon et all [2010l 
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2012t iTilvi eralll2010t IKrue et al.ll2012D. or fuh SEP 
fitting to broadband photometry (e.g.,lFinkelstein et al.l 
[20T0ll20Tlt IMcLureet aniMoL mm have required a 
large investment of observing time on HST and/or large 
(> 4-m) ground-based telescopes and have detected only 
a handful of candidates. 

Gravitational lensing by galaxy clusters aids in the 
detection of these sources by magnifying the back- 
ground source population. A number of surveys toward 
known lensing clusters ha .ve revealed a few high-redshift 
galaxy candidates (e.g., iKneib et al.l l2004t iPello et all 
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ing makes background sources easier to detect by map- 
ping regions from the source plane to a larger area in the 
lens plane, there is a corresponding loss of volume probed 
in the source plane. These two effects (gain in depth, 
loss in volume) affect the efficiency of searches for high- 
redshift galaxies in different ways. Addressing the ef- 
fects of lensing magni fication on the detectability of high- 
redsh ift sources (e.g., IBouwens et al.l [20091 : iMaizv et al.l 
12010( 1 is subject to uncertainties in the source properties, 
such as the poorly-constrained shape of the z > 7 galaxy 
luminosity function (LF), the lens properties, and the 
observational strategy. The efficiency of massive cluster- 
scale halos at magnifying a background source popu- 
lation is affected by halo mass, ellipticity, orientation, 
and substructure, as well as b y uncorrelated structures 
along the line of s i ght (LOS: IWambsganss et al.l [20051 : 
IHilbert et al1[2007t IPuchwein fc Hilbertll2009D . 

The lensing properties of a given LOS, or "beam", de- 
pend not only on the total integrated mass along the 
beam, but also on the distribution of that mass. Ideally, 
that mass should be distributed in a way that maximizes 
the etendue, the area in the source plane over which sig- 
nificant magnification occurs. Selection of cluster-scale 
lenses through X-ray observations and weak lensing stud- 
ies, among other methods, has often yielded fields where 
a single massive structure dominates the beam, while 
any smaller LOS structures are ignored. Yet, beams 
could exist where multiple cluster-scale structures max- 
imize the contribution to the lens potential, enhancing 
the detectability of lensed background sources over that 
of a single lensing cluster, even if the total integrated 
mass in the beam is the same. 
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A few studies have explored the theoretical lensing 
prope rties of lines of sigh t containing multiple clus - 
ters HCrawford et all \T9m iBertin fc Lombardil [200l . 
often focusing on cluster-cluster strong lensing (e.g. 
ICoorav et al.l Il999f) . but few such configurations have 
been observed and characterized in terms of the prop- 
erties that maximize their lensing strength. Possi- 
ble chance al ignments in cluster o bserv ations were also 
sugges ted by iMolinari et al.l (|1996[) and iWang fc Ulmeii 
(|1997| ). Serendipitous discoveries have since shown that 
chanc e alignmen ts of cl usters exist in the observable uni- 
verse ()Blakesleell200lt lAthreva et all 120021: [zTtrin et al.l 
120121) . However, the lensing properties of such multi- 
cluster beams relative to single clusters have not been 
well-studied, and the likelihood of these alignments is 
not well-known. Distributing the mass in a beam among 
multiple projected lensing halos may lead to beneficial in- 
teractions of the lens potentials and thus higher etendue. 
If these beams are better at magnifying distant sources 
than typical lensing clusters, it is important to identify 
the characteristics that make them powerful lenses and 
to quantify their frequency in the universe. 

The physical properties of the halos in multi-halo 
beams can also impact their ability to lens high-redshift 
sources. Several past studies have examined the effects 
of halo concentration, shape, and substructure, although 
in the context of lensin g z ^ 1 — 2 galaxi es into arcs. 
IDalal et al.l (|200l and IMandelbaum et al.l |2009) find 
a strong correlation between halo mass and the lens- 
ing cross section, consistent with arcs in galaxy clus- 
ters and not in lower mass systems. As a result of 
numerical studies of halo concentration, ellipticity, and 
triaxiality, several groups (e.g., [B artclm aiin et al.l 119951 : 
IMeneghetti et al.ll2003HOguriet all 1200^ find that the 
lensing cross section increases as the mass distribution 
of the halo becomes more centrally concentrated (either 
through a steepening of the inner density slope or an in- 
crease in the projected halo concentration). These stud- 
ies have focused on single halos and generally do not seek 
to identify the properties of the very best and realistic 
lensing configurations for very high-redshift galaxies. 

In this paper, we perform a Monte Carlo analysis to 
model beams containing multiple cluster-scale halos and 
evaluate their lensing properties in comparison with sin- 
gle cluster lenses of equivalent total mass. By controlling 
for total mass, we isolate the gain due purely to the mass 
configuration — the number of halos and their physical 
properties. We apply physically-motivated priors to en- 
sure that our beams are realistic. 

We describe the details of our lensing analysis and out- 
line our methodology for characterizing the best lensing 
beams in § [21 In § [31 we describe the models we use 
to represent various mass configurations and test their 
lensing properties. We present our results in § [H We 
discuss the interpretation of our results in § [5l as well as 
possible caveats and sources of uncertainty. In addition, 
we discuss the implications of our findings on the detec- 
tion of high-redshift sources and the likelihood of find- 
ing such mass configurations in the universe. Our main 
conclusions are summarized in § [Bl Throughout this pa- 
per, we assume a ACDM cosmology with £7™ = 0.274, 
Qb = 0.045, Ha = 0.726, and Hq = 71 km s'^ Mpc"!. 



2. CHARACTERIZING THE BEST MASS 
CONFIGURATIONS FOR LENSING 

The lensing properties of a particular mass configura- 
tion for a high-redshift source plane can be infiuenced 
by several physical characteristics. The total mass along 
the line of sight clearly is important, as well as the man- 
ner in which that mass is distributed, both in projection 
and redshift. In addition, the structural properties of the 
halos containing the mass will have an effect on the lens- 
ing. In this section, we describe our lensing calculations 
and our metric for evaluating the lensing properties of 
our simulated beams, as well as introduce our tests for 
looking at the effects of the various physical properties 
of their mass configurations. 

2.1. Lensing Analysis 

All of our lensing calculations are performe d using 
an upd ated version of the software described in iKeetonI 
(|2001aD . Here, we describe methodology for determining 
the lensing properties of a given mass configuration. We 
select a source redshift of z^ = 10, which is close to the 
current limits of what is currently fe asible with existing 
observatories (jBouwens et al.l [20 Hal) . 

To study the lensing properties of various mass config- 
urations on a z = 10 source plane, we construct a uniform 
5'x5' grid of small circular sources at 1" intervals on the 
source plane. In principle, we could use elliptical sources 
()Keeton|[2001b[) . but this adds an additional complica- 
tion, so we assume circular sources. We then trace these 
sources through the mass configuration. The result is a 
distribution of image positions that will always contain 
a number of points greater than or equal to the origi- 
nal source plane grid due to multiple imaging of sources 
within the lensing caustics. For each source, we tabulate 
the images associated with it, as well as the elements of 
the magnification tensor, M(a;), for each image. The to- 
tal image magnification at x is just the determinant of 
M{x). 

The axis ratio of the image (assuming an intrinsi- 
cally circular source) can be approximated as the ra- 
tio of the eigeiivalues of the magnificat ion tensor (e.g., 
ISchneideret all [1991 IFedehet al.l [2001 . a/b w A1/A2, 
where a and b are the major and minor axes of the 
lensed image and Ai and A2 are the maximum and min- 
imum eigenvalues of M(a;), respectively. This approxi- 
mation should be valid for small circular sources, but can 
break down for extended sources when the magnification 
tensor changes appreciably over the area of the source 
(Fedcli ct al. 2006). 

The magnification tensor at each image position is cal- 
culated numerically, accounting for the effects of multiple 
lensing planes. We provide a brief discussion of the lens- 
ing formalism behind these calculations since we are of- 
ten dealing with multiple halo s at di fferent redshifts (see 
also lKocha-nek fc ApostolakisI (|1988D and McCuUy et al. 
in preparation). Angular positions in the source plane 
and image plane are denoted u and x, respectively. 

For simplicity, first consider the case of two halos at 
the same redshift so they can be treated as being in the 
same lens plane. The lens potential is the superposition 
of the potentials of the two halos. 



(t>totix) = (j)i{x)+(l)2ix). 



(1) 
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The lens equation is 

U = X - V^tot 

= X — ai{x) — a.2(x), 



(2) 



where cti = V(/)i and 0.2 = V02 are the deflections for 
the two halos. The inverse magnification tensor is 



M(a;)-^ = 



1 — Kl 



- K2 



-Id 
Isl 



7c2 



1 — Kl 



-7s 1 - 7s2 
- K2 + 7cl + 7c2 
. . (3) ■ 

where k indicates the convergence (the isotropic magnifi- 
cation term) while 7^ and 7^ indicate the two components 
of shear (the anisotropic magnification term). In general, 
these three quantities will have similar amplitudes — not 
necessarily identical, but within a factor of a few of each 
other. The magnification at a given image position is 
then 



/Li = det(M(a;)) 
= [(1 



1 



det(M(a;)-i) 



Kl - K2 ) 



(7c 



-7c2j 



(7. 



■Is2) 



(4) 



The convergence and shear terms from both halos enter 
the equation at the same lowest order. 

Consider a scenario where the two halos are separated 
enough that their caustics are distinct. We would like to 
examine how the lens potentials of the halos interact with 
each other and affect the magnification in regions near 
one of them. Suppose we are reasonably close to halo #1. 
Equation|4]indicates that the change in the magnification 
across the field due to the presence of halo #2 can be due 
to both its density profile and shear effects. If halo #2 
is far enough away that its convergence and shear are 
relatively small, we can do a Taylor series expansion in 
the quantities from halo #2, 



^ « ^1 + 2ll\ [(1 - Ki)k2 + 7cl7c2 + 7sl7s2] + 



(5) 



-ll. 



7^2] ^ is the magnification 



where ^1 = [(1 — Ki)^ 

from halo #1 alone. The extra term in Equation ([5]) 
does not have the form of the magnification expected 
from halo #2 alone. This is the statement that — unlike 
the lens potentials, defiections, convergence terms, and 
shear terms — magnifications do not simply add linearly 
because magnification is a nonlinear function of the po- 
tential. This term is not necessarily positive since 7c and 
7s can have either sign, although in practice the presence 
of halo ^2 generally increases the magnification rather 
than decreases it. 

Now consider the case of two halos at different red- 
shifts. Let halo #1 be in the foreground and halo #2 be 
in the background. Now the lens equation has the form 



u = X — a.i{x) — 0:2(2; — /3i2Q;i(a;)), 



(6) 



where /3i2 is a ratio of angular diameter distances be- 
tween combinations of the observer, the two lens planes, 
and the source plane: 



/3i 



D12DS 
D2D1, 



(7) 



Note in Equation |6] that a.i appears inside the argument 
of a.2- This occurs because the location where a light 



ray intersects plane #2 is shifted from the position of 
the light ray in the image plane due to the presence of 
the foreground lens. The inverse magnification tensor is 



m{x)-^ = 



1 - (ki +7ci) 
-7s 1 1 



K2 + 7c2 7s2 
7s2 K2 - 7c2 



-7si 
(ki - 7ci) 

- -/5i2(ki +7cl) 

-Pl2lsl 1 



-/3i27si 
/3i2(ki -7ci) 

(8) 



This is actually even more complicated than it appears 
because (like 0L2) the functions K2, 7c2, and 7s2 are all 
evaluated at the position x — f3i2a.i{x). The expression 
for /i is similarly complicated. 

The formali sm can be extend[ed to an arbitrary numbe r 
of lens planes ()Schneider et al.lll992t iPetters et aDl200ir ) . 
The inverse magnification tensor for the jih. lens plane 
of a given line of sight can be calculated through the 
recursion relation 



i-i 



m{x)f = \-y^p,,TM{x)-;\ 



1=1 



(9) 



where I is the 2x2 identity matrix, M(a;)j^ ^ = I, and 
Ti is the "shear tensor" for plane i, 



Is'i. Hi - 7c. 



(10) 



The inverse magnification tensor for the entire mass dis- 
tribution is given by Equation [S] evaluated at the source 
plane, M(a;)j^, and the image magnification and axis ra- 
tio can then be evaluated. We will revisit this formalism 
when interpreting our results. 

2.2. Quantifying the Lensing Properties of Different 
Mass Configurations 

There are a number of quantities we can use to deter- 
mine how good a particular line of sight is for lensing. 
The relative merits of each quantity depends on our ul- 
timate objective. For the purposes of this work, we are 
primarily interested in finding the lines of sight that are 
most effective at magnifying very high-redshift sources 
(z ^ 10) into detectability. 

A simple and useful metric for evaluating the lens- 
ing power of a particular mass configuration is the cross 
section for magnifying a source above some minimum 
threshold value. This cross section, which we denote 
(jf^, is the area in the source plane over which a lensed 
source's brightest image is magnified above the chosen 
threshold magnification, a/^^ is a convenient quantity be- 
cause it makes no assumptions about the properties of 
the source population (aside from its redshift) or the 
observational setup, but is dependent only on the mass 
configuration. We choose to define af^ in terms of source 
plane area rather than image plane area because the for- 
mer is a better metric for evaluating the quality of a par- 
ticular line of sight in terms of maximizing the potential 
number of high-redshift detections. Although area in the 
image plane with high fj, may be a relevant quantity, its 
usefulness depends on the specific details of the observa- 
tions designed to study these lines of sight due to source 
size effects, and is therefore less general. 
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We define tr^ in terms of tlie magnification of the 
brightest source image rather than total magnification 
because it is the detectabihty of individual images that 
interests us. If a source is multiply imaged such that 
it has a high total magnification but no individual im- 
age is magnified enough to be detectable, it is of no use. 
Throughout this paper, we define tr^ as the source plane 
area where a source's brightest image has a magnification 
/i > 3. This value is somewhat arbitrary, although there 
are indications that the number of high-redshift detec- 
tions in lensing fields, while relatively fiat as a function 
of magnification, broadly peaks in regions with /i ~ 3 
(Amnions et al. in preparation). We ultimately find that 
our general conclusions are not affected by our choice of 
magnification threshold, and we discuss the implications 
of changing this threshold in Appendix |^ Formally, cr^ 
is defined as 



(J^= n (max(^) - ^f) 
Jn 



dfl^ 



(11) 



where f2 is area in the source plane, H is the Heaviside 
step function, max(/i) is the highest magnification of any 
individual image associated with a given source position, 
and /it is the threshold magnification, which we choose 
to be fj.t = 3. a^ can be thought of as etendue — the 
area in the source plane that would get magnified above 
a certain value of interest. 

Figure [1] shows a^ as a function of mas s for a single 
spherical NFW halo (jNavarro et al.lll996D at z = 0.5, 
with structural properties at each mass given by the me- 
dian relations of I Zhao et al.l (|2OO90 . The halo mass is de- 
fined as M200, the mass enclosed within a sphere whose 
mean density is 200 times the mean matter density at 
that redshift. The mass range is chosen to represent re- 
alistic sizes for massive cluster lenses. As expected, the 
cross section increases monotonically with mass, chang- 
ing by >3 X for a corresponding factor of 3 change in mass 
over this range. This refiects the fact that more massive 
halos are generally better lenses. This figure also gives 
a sense of the value of a^ expected for massive clusters 
and how it changes with the addition of more mass. We 
later use this model to roughly compare gains in cr^ aris- 
ing from the mass configuration alone to those from the 
addition of more total mass. 

3. LINE-OF-SIGHT MASS MODELS 

The halos used in our analysis are parameterized as 
NFW profiles, an assumption that facilitates comparison 
with past and on-going work on massive cluster lenses 
and is supp orted by se veral observati o nal studies (e.g., 
Kneibet al.l [20 03: Kat gert et al.l[200l iComerford et alJ 
20061: iShu etar.li2008iV We neglect the effects of substruc- 
ture and baryonic matte r in our calculations, although 
we comment on this in S 15.21 



3.1. Priors on Halo Properties 

We generate halos with redshifts drawn from a distri- 
bution in the range 0.1 < z < 0.9. We choose this range 
because we expect the majority of massive cluster-scale 
lenses to be at these low-to-intermediate redshifts due 
to the evolution of the halo mass function. At very low 
redshift (z < 0.1), there is little volume and the sur- 
face density of cluster halos becomes small (assuming 
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Fig. 1. — The lensing cross section, cr^, as a function of mass for a 
spherical halo at z = 0.5. The cross section increases monotonically 
with mass, reflecting the fact that more massive halos are better 
lenses. 



a shallower-than- isothermal mass profile), making these 
halos less effective lenses. At higher redshifts {z > 1), 
the lensing geometry works against gains in magnifica- 
tion because the angular diameter distance ratios among 
the observer, lens, and source become unfavorable for 
lensing high-redshift sources. Furthermore, cluster-scale 
halos are more rare and the evolution of cluster prop- 
erties has not been well constrained by observations at 
these redshifts. Our redshift range encompasses the the- 
oretical peak in the lens redshift distributio n for clus- 
ters l ensing giant arcs {z ~ 0.3 — 0.4; Bartcl mann et all 
Il998f ). as well as the observed peak in the redshifts of 
arc-produc ing lensing clusters f r om the RCS2 datase t 
(z > 0.6: IGladders et al.l [200I IHennawi et~all [200I . 
The redshift distribution from which we draw is weighted 
by the number of halos with masses greater than 10^^ Mq 
expected within a comoving volume element defined by 
the solid angle of our fiducial aperture of 3.'5 (see below). 
This procedure ends up populating the full redshift range 
(0.1 < z < 0.9), although not uniformly. Our procedure 
includes the possibility that multiple halos lie at the same 
redshift, as would subclusters in clusters of galaxies, al- 
though we do not account for physical interactions among 
the halos as might occur in subcluster mergers. 

The number density of halos at a given mass and red- 
shift is determined using the pub l icly-a vailable halo mass 
fun ction code of 'Tinker et al.l (|2008D . which assumes 
the lEisenstein fc Hu (1998) matter transfer function as 
the default. In this calculation, we also assume our 
fiducial cosmological parameters, along with erg = 0.8, 
Scrit — 1-686, and a primordial power spectrum index of 
0.963. We define halos to have an overdensity of A = 200 
relative to the mean background density at that redshift. 
We assign each halo a mass drawn from the mass function 
at its given redshift, requiring that the halo has a mass 
of at least IO^'^Mq, roughly the mass of a Virgo-scale 
galaxy cluster. 

We perform a Monte Carlo analysis to generate lines 
of sight that contain N halos, where N is the number of 
distinct halos that we are testing in that particular set 
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of trials. We examine a control sample of beams con- 
taining A^ = 1 halo (§ 13. 2p . then explore beams with 
multiple halos for comparison fij 13.31 and ii l3.4p . For each 
trial, we randomly draw redshifts for N halos using the 
procedure described above. We then randomly assign 
masses to these halos using the mass function at each 
halo redshift until we find a realization where the to- 
tal mass across all N halos is within 1% of some fidu- 
cial beam mass. We choose this fiducial mass to be 
2 X 10^^ Mq which is roughly the virial mass of massive 
clusters like Coma (Kubo et al. 2007) and AbeU 1689 
(|Umetsu fc Broadhurstll2008f ). 

The angular positions of the halos in a particular re- 
alization are drawn randomly to approximate the dis- 
tribution of physically uncorrelated halos in such a line 
of sight within a fixed aperture. We choose an aperture 
size of 3.'5 in diameter, which is a separation large enough 
that halos with our chosen range of properties generally 
behave as independent lenses and do not gain from in- 
teractions with other halos (see § 14. ip . We place one 
halo at the origin and randomly distribute the remain- 
ing A^ — 1 halos uniformly across a disk of radius 3'. 5. 
We then calculate the centroid and diameter of the mini- 
mum enclosing circle of the distribution. If the diameter 
is less than 3'.5, we retain that configuration, otherwise 
we repeat the procedure. After finding a suitable config- 
uration for a particular trial, the field across which we 
perform our lensing analysis is re-centered on the cen- 
troid of the minimum enclosing circle. 

For beams containing multiple halos, this procedure 
assumes that the positions of cluster-scale halos on the 
sky are uncorrelated on scales smaller than a few arcmin- 
utes. In detail, there may be correlations between halos 
at small redshift separation due to large scale structure 
(e.g. filaments, superclusters). However, accounting for 
this would require detailed knowledge of the evolution 
of the cluster correlation function on ~arcminute scales 
with redshift, which is not well-constrained. Such corre- 
lations will be investigated in a future work using cos- 
mological simulations (French et al. in preparation), but 
we ignore this effect for our analysis. 

To determine the relevant physical properties of the 
halos for our lensing calculations, we use the publicly- 
available h a lo ev olution code based on the models of 
I Zhao et all ()2009f ). Using this code, we determine the 



median virial radius 



concentration parameter Cy 



scale radius r^, and characteristic density ps as a func- 
tion of halo mass and redshift and apply these to our 
halos, assuming the same cosmological parameters as de- 
scribed above. We note that recent results suggest a 
possible reversal in the halo mass-concentration relation 
at th e cluster scale (e.g., Klypin et al. 2011; Prad a et al.l 
1201 It ), which we do not account for. However, such a 
trend would affect both our control sample and multi- 
halo beams in similar ways, so the relative properties of 
these lines of sight should not be significantly affected. 
We include an intri nsic scatter in Cyir of CTiogc = 0.14 
(jBullock et al.l [20011 ). which affects the values of r^ and 
Kg (the characteristic surface density) that go into the 
lensing calculation but keeps virial mass and radius fixed. 
Lensing clusters have been found to have an orienta- 
tion bias such that the major axis of lensing-selected 
halos tends to be preferentially aligned with the line 
of sight from the observer. This has the effect of 



increasing t heir project ed concen tration in the plane 
of the sky (iCorless fc K ing 200 71: iHennawi et all 120071: 
iBroadhurst fc Barkanal [2008t lOguri fc BlandfordI I2009D . 
We are not looking specifically at lines of sight that have 
been predetermined to be good lenses, but triaxiality and 
projection effects can certainly affect the lensing prop- 
erties of these mass configurations. For this reason, we 
assume that our halos are triaxial, with axis ratios drawn 
from Gaussia n distributions det ermined from simulations 
performed bv lShaw et al.l (|2006l ). We use the results from 
their simulations that include mass structures just be- 
yond the virial radii of their halos, which give axis ratios 
oib/a = 0.79± 0.121 and c/a = 0.678± 0.113, where a, &, 
and c are the major, intermediate, and minor semi-axes 
of the halo. We then give our halos a random orien- 
tation in a spherical coordinate system and project the 
mass distributions onto the plane of the sky to get their 
projected ellipticities and position angles. We can also 
determine the effects of this projection on the relevant 
physical quantities for lensing (see Appendix |B] for de- 
tails). We note that our assumed NFW halos have ellip- 
tically symmetric density profiles (as opposed to ellipti- 
cally symmetric lens potentials) for which the potentials 
do not have analytic solutions and are instead computed 
numerically. 

Any values of lensing quantities given here are spe- 
cific to the mass models with the properties used in our 
analysis. In general, the diversity of possible mass con- 
figurations in the universe likely results in a large range 
of possible lensing properties and implied effects on the 
detection of high-redshift sources. We have imposed sev- 
eral physically-motivated priors to reflect realistic mass 
distributions, but we caution against using the absolute 
numbers derived in this analysis, instead focusing on 
comparing the relative values of the lensing properties 
of the multiple-halo lensing beams to that of the single- 
halo case. 

3.2. Single Halo Control Sample 

We generate a sample of fields containing a single halo 
as a control sample with which to compare the multi-halo 
mass distributions. These fields contain a single NFW 
halo with a virial mass within 1% of our fiducial mass of 
2 X lO^^Af0 with properties described in § 13.11 A source- 
plane magnification map for a typical one of our Monte 
Carlo trials is shown in the top left panel of Figure [21 
with the corresponding image-plane map shown in the 
top left panel of Figure [3l Cross-sections are evaluated in 
the source plane. We will also investigate the properties 
of the best mass configurations that are ranked within 
the top 10% by a^, so we plot analogous maps in the top 
left panels of Figures [4| and [5] for one such configuration. 

3.3. Two Halos 

Our first multi-halo test is to look at lines of sight with 
a total mass of 2 x IO^^Mq divided between two separate 
cluster-scale halos. Sample magnification maps for one 
of our typical trials are shown in the top right panels 
of Figures [2| and [H with the same plots for a top 10% 
trial show in the top right panels of Figures [D and [5j 
Although the total mass in the beam remains fixed to 
within 1%, we can quantify the differences in their lensing 
properties compared to the single-halo beams of similar 
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containing a single halo (top left), two halos (top right), three halos (bottom left), and five halos (bottom right). The colors correspond 
to the magnification of the single most magnified image of a small source at each position on the map. The color map is logarithmic in 
the magnification. The values of (J^^ (arcsec ) for these particular configurations are indicated in the top right corner of each panel. These 
median configurations show that increasing the number of halos results in an increase in Cfi . 



total mass, as well as investigate the effects of varying 
other properties such as the halos' angular separation, 
relative mass ratio, and redshift distribution. 

Although we account for the halo mass function and its 
redshift evolution when assigning relative halo properties 
in our analysis, we are, by construction, forcing lines of 
sight to have both a fixed number of halos and a fixed 
amount of total mass, not accounting for the likelihood of 
finding a high-mass bea m containing multiple halos pro- 
jected nearby on the sky. lOguri fc Blandfordl (|2009P ) find 
that chance superpositions of two cluster halos that to- 
gether create a large Einstein radius are rare compared to 
single clusters with similar lensing properties. However, 
we are by design investigating the properties of lines of 



sight that contain rare mass configurations that are fa- 
vorable to lensing. We address the likelihood of these 
chance projections existing in the universe in § 15.41 The 
frequency of such fields will be analyzed in more detail 
using cosmological simulations in a future work (French 
et al. in preparation). We note, however, that we have 
identified several such lines of sight in the SDSS but will 
leave analyses of their properties to future work (Am- 
mons et al. in preparation). 

3.4. More Than Two Halos 



We generalize the tests performed in § 13.31 to allow for 
the total mass in the beam to be redistributed among 
three or more halos projected along the LOS. For a trial 
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with TV distinct halos, we draw N halos until we find a 
beam with a total integrated mass within 1% of our fidu- 
cial single halo. We test lines of sight containing A^ = 3 
and N = 5 halos. Sample magnification maps for typical 
3-halo and 5-halo cases are shown in the bottom row of 
Figures[5]and[31 with a top 10% case shown in the bottom 
row of Figures |4] and [S] Further trials of beams contain- 
ing even more distinct halos are possible, but calculating 
their lensing properties becomes computationally inten- 
sive. Since these lines of sight are increasingly rare in 
the real universe as more cluster-scale halos are added, 
we do not explore beams with more than five halos. 

4. RESULTS 

We present the results of our analysis of the lensing 
properties of the single-halo and multi-halo lines of sight 



and look into the physical properties driving the differ- 
ences between them. In addition, we assess the frequency 
of various properties of the mass configuration and inves- 
tigate the parameters that lead to the very best lensing 
configurations, independent of frequency. 

4.1. The Effect of Multiple Halos On Lensing 

The results of our Monte Carlo analysis are presented 
in Figure |6l where we show the distribution of a^ for 
all 5000 trials of each LOS mass distribution. These re- 
sults indicate that the multi-halo beams are, on aver- 
age, better at magnifying high-redshift galaxies into de- 
tectability than single-halo beams of similar integrated 
mass. Furthermore, this effect becomes more prominent 
as the mass is redistributed over more distinct halos. The 
statistics of the distributions are quantified in Table [TJ 
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Fig. 4. — Sample z = 10 source plane magnification maps over a 300" X 300" field for lines of sight ranked near the top decile in cr^ 
with a single halo (top left), two halos (top right), three halos (bottom left), and five halos (bottom right). The colors correspond to 
the magnification of the single most magnified image of a point source at each position on the map. The color map is logarithmic in 
the magnification. The values of ct^ (arcsec^) for these particular configurations are indicated in the top right corner of each panel. The 
increase in Cfi as the number of halos increases is similar to the result shown in Figure [2] 



The typical improvement in the average lensing cross 
section of the 2-halo beams over the single-halo beams 
of equivalent mass is roughly ~ 15%. The improvement 
for the 3-halo and 5-halo beams are roughly ~ 40% and 
^ 100%, respectively. These gains are specific to our 
chosen fiducial mass and parameter ranges and will differ 
for various mass and separation range, but the overall 
eflfect remains similar. 

It may be the case that splitting a single massive halo 
into multiple smaller halos increases the lensing cross sec- 
tion simply by virtue of covering a larger area in the sky. 
It is not immediately obvious that it is the combined ef- 
fect of multiple halos nearby in projection that is creating 
this effect. To check this, we run a test where for each 2- 



halo trial, the two halos are analyzed as if they were each 
individual clusters, and their total lensing cross sections 
are summed. Figure [7] shows histograms of tr^ for these 
trials compared to both the single-halo beams and the 
2-halo beams. Table [T] also shows the mean and median 
values of the cross section for these trials. The results 
show that simply dividing a single halo into two smaller 
independent halos does not result in the same gain that is 
seen when it is split into two halos acting in combination 
along the same line of sight. In fact, doing so actually 
decreases the lensing power of the mass distribution. It 
is the combination of the two halos in projection acting 
together that results in the increase seen in Figure ID We 
discuss this eflFect further in S 15.11 
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Number of Halos 



TABLE 1 
Lensing Properties 



{a-fj,)'^ (arcsec ) 



Median(a'p) (arcsec ) 



Top 10% Decile in cr^ (arcsec ) 



1 
2 
3 

5 

2 (individual)*^ 



1965 ± 5 
2268 ± 6 
2747 ± 8 
3939 ±11 
1835 ± 5 



1990 
2262 
2708 
3859 
1868 



2300 
2799 
3470 
4972 
2190 



''Errors given are standard errors on the mean. 

''Each of the beams in the 2-halo trials were split into individual halos and the lensing cross sections of each were summed. 
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Fig. 6. — Histogram showing comparisons of a^ for the single 
halo lines of sight to those of the multi-halo lines of sight from 
our Monte Carlo analysis. 5000 beams of each configuration were 
generated for each case. The beams containing multiple halos have, 
on average, larger cross sections than the single-halo beams. 
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Fig. 7. — Same as Figure [6] but comparing only the single-halo 
and 2-halo lines of sight, as well as the sum of the quantities in 
each trial determined from each of the individual halos separately 
in the 2-halo cases. The beams containing two halos have, on 
average, larger cross sections compared to the two halos considered 
separately. 

Additionally, one can compare the gain in cross sec- 
tion due to configuration to the gain due to an increase 
in total mass, as shown in Figure [TJ While this basic 
comparison only represents the scaling for one simplified 
scenario, a general sense of how much the configuration 
can afi^ect the cross section compared to a change in the 
total integrated mass can be inferred. For example, com- 
paring the mean cross section in the single halo case to 
that of the 5-halo case shows an increase of a factor of 
~2, which is roughly comparable to changing the total 
mass of the single halo by a factor of ^-^2. This demon- 
strates that the effect of configuration is non-negligible 
compared to the effect of total mass, and may even be 
comparable over this range in total mass. 



4.2. Halo Properties Important for Lensing 

It is important to determine what parameters drive 
correlations between the physical properties of the multi- 
halo mass distributions and their lensing properties. We 
focus on the 2-halo beams for this discussion, as many 
of the relevant physical properties of these lines of sight 
can be represented by a single quantity, and 2-halo beams 
will be more common than lines of sight containing three 
or more cluster-scale halos. 

The results of § 14.11 show that the lines of sight com- 
prising multiple halos tend to be better lenses than the 
total of each of their constituent halos considered indi- 
vidually, assuming fixed total mass. This suggests that 
for a given set of halos, there is an optimal configuration 
in the plane of the sky that maximizes the lensing signal. 
For the 2-halo case, this corresponds to an optimal an- 
gular offset given two particular halos. In the left panel 
of Figure [H we plot (jf^ for each trial as a function of the 
angular separation between the two halos. Due to the an- 
gular positions of the halos relative to one another, there 
is a higher overall frequency of beams at larger projected 
separations, as we would expect if the halo positions on 
the sky are uncorrelated. Despite the fact that this vi- 
sualization marginalizes over all other parameters, it is 
clear that cr^ rises gradually, peaks at a separation of 
^100", and then decreases at higher separations. 

The dropoff at large separation appears to be sensi- 
ble in the context of the results of Figure U\ since we 
are approaching the limit where the regions around each 
halo can be treated as independent fields. The peak at 
nonzero separation is more interesting, and can be ex- 
plained by the distribution of surface mass density when 
multiple mass components are spread apart. We discuss 
the interpretation and implications of this result further 
in §0 

We perform a similar test to examine the effect of vary- 
ing the relative masses of the halos in the beam. The 
right panel of Figure |8] shows cr^ as a function of the log 
of the ratio of the two halo masses (defined such that the 
mass ratio is always < 1). The trend shows that mass 
ratios closer to unity results in larger cross sections. As 
the mass ratio becomes more extreme, we are approach- 
ing the limit of the single-halo case, so this result makes 
sense intuitively. As an additional check of these results, 
we run 2-halo trials where we both increase and decrease 
the minimum mass allowed for our halos (see S I5.2|) . Both 
of these tests show the same trend of tr^ with mass ratio, 
and the mean cr^ across all trials increases as the min- 
imum mass (and thus, the mean mass ratio) increases, 
supporting our conclusion. 

We note that Figure |8] indicates a large number of trials 
with extreme mass ratios where there is a single domi- 
nant halo. Due to the steepness of the halo mass function 
above IO^^Mq within our redshift range, drawing N halo 
masses that add up to our fiducial mass is much more 
likely to select one massive halo along with A^ — 1 halos 
close to our minimum mass of IO^^Mq than it is to select 
A^ similarly massive halos, which results in this skewed 
distribution. One may wonder in this case whether the 
lensing properties of the multi-halo beams are dominated 
by the physical properties of the more massive halo. To 
check this, we look at how the properties of the single- 
halo beams affect the lensing cross section and compare 
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Fig. 8. — cr^ as a function of angular separation (left) and mass ratio (right) of the hales in the 2- halo trials. The grey points represent 
individual trials, and the black points are averages in each bin. The error bars represent standard errors on the mean. Left: a^ rises 
with separation and peaks around ~100", then drops off at larger separations. This is consistent with the understanding that at large 
separations, the halos approach the limit where the regions around them can be treated as independent lensing fields. Bight: The mass 
ratio is defined as the mass of the smaller halo divided by the mass of the larger halo. Mass ratios closer to unity produce larger lensing 
cross sections, with cr^ generally decreasing with smaller mass ratio. At low mass ratio, we approach the limit of the single-halo fields, 
which have lower cross sections on average than the 2-halo case. 



any trends or features to those in the multi-halo beams 
where only the properties of the more massive halo are 
considered. 

In the top row of Figure IH we plot cr^ for the single- 
halo lines of sight as a function of various fundamental 
parameters. We see that halos at redshifts higher than 
z ^ 0.3 tend to be better lenses, which is due to the 
fact that the critical density for lensing, Ec, has a broad 
minimum centered at z ~ 0.8, resulting in an optimal 
geometry for lensing z = 10 sources. We also see that 
halos with higher projected concentration, c', tend to 
have larger cross sections, which is expected due to the 
increase in surface density that this implies, c' is deter- 
mined both by the intrinsic concentration of the halo, 
as well as density enhancements due to triaxiality and 
projection effects (see Appendix iBt . This results in ad- 
ditional scatter in c', with halos that have their major 
axis aligned with the line of sight tending to have larger 
concentrations, and therefore larger cross sections. We 
use uj to denote this angle, which is defined to be zero 
when there is perfect alignment between the major axis 
and line of sight direction. The effect is small in compar- 
ison to the intrinsic concentration of the halo, but we see 
in the third panel of Figure |9] that it is not negligible. 

Other physical parameters do not show a statistically 
significant correlation with afj,, so we do not plot these 
quantities. These parameters include projected elliptic- 
ity and triaxiality, as well as relative ellipticity position 
angles between both halos or between the more massive 
halo and the angle to its companion. Given that there is 
a trend with w, one might have expected a possible trend 
with halo triaxiality. However, because w is defined by 
the direction of the major axis, regardless of halo shape, 
this is not the case. Likewise, multiple combinations of 



axis ratios and viewing angles can result in similar pro- 
jected ellipticities, so there are no obvious dependencies 
on ellipticity alone. 

We also plot cr^i as a function of the projected scaled 
surface density for NFW halos, k'^ (see Appendix IBJ) . in 
the upper rightmost panel of Figure [9] k'^ is dependent 
on other fundamental parameters, but we show it here 
because this relation shows an interesting shape. Very 
low values of k'^ correspond to small cr^, which occurs 
when the halo is at a low redshift such that it has a large 
angular size and therefore a small surface density. One 
might expect that the cross section would monotonically 
increase as the surface density increases, but instead, cr^ 
reaches a point where it plateaus and eventually turns 
over. This unexpected behavior has to do with our def- 
inition of the cross section. As k'^ increases, the area 
in the source plane interior to the inner caustic grows. 
Any sources that fall within the caustic get multiply im- 
aged, but the maximum magnification of any individual 
image is sometimes lower than our threshold magnifica- 
tion, /it. Therefore, this area does not count toward the 
cross section as we have defined it. Obviously, different 
magnification thresholds will change the turnover point, 
but adding more surface density in general will increase 
this area but will not increase a^ once a certain point 
has been reached. We note that the traditional strong 
lensing cross section, defined as the area in the source 
plane where a source is multiply-imaged, does increase 
monotonically with n'^ as expected. 

We examine whether these results also hold true for 
the multi-halo distributions where only the properties of 
the more massive halo are considered. In the bottom row 
of Figure El we plot cr^ for the 2-halo lines of sight as a 
function of the same parameters for the more massive 
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Fig. 9. — o-^ as a function of (from left to right) redshift, projected concentration, angle between the major axis and the line of sight 
(denoted uj), and projected scaled NFW surface density k^ of the halos in the single-halo sample (top) and the more massive halo in each 
trial in the 2-halo sample (bottom). Other parameters, including projected ellipticity and triaxiality, are not plotted because they do not 
show a statistically significant trend, o-fj, is larger in general at redshifts z > 0.3 compared to lower redshifts, but is fairly constant with 
redshift on average beyond that point, cr^ also increases with larger projected halo concentration and decreasing uj. Cfi rises sharply with 
K^ at low values, then plateaus due to the way we have defined the cross section in our analysis (see text for details). The results in the 
two-halo case show that the properties of the more massive halo affect the cross section in the same way as the single-halo case because 
most trials have a single dominant halo. We note that c' and re^ are the projected concentration and scaled surface density after accounting 
for triaxiality and projection effects (see Appendix IB)| . 



halo in each triaL The properties of the dominant halo 
affect the cross section in a similar way as in the single- 
halo case, as we would expect when one halo dominates 
the mass distribution. While 1:1 mass ratios generally 
produce higher ct^, even the far more numerous > 10:1 
mass ratios can generate high ct^ . One implication of Fig- 
ure E] is that in studies of known strong lensing clusters 
where LOS structures are typically ignored, any mas- 
sive structures in the beam that go undetected are more 
likely to be small clusters rather than one with a mass 
equivalent to the main lens. As such, not accounting for 
these structures is unlikely to result in a large error on 
the cross section of these beams. However, the effects of 
these structures on the derived 2-D magnification maps, 
which we do not analyze here, may result in systematic 
uncertainties that will cause the magnification errors in 
the field to be underestimated. 

4.3. Properties of the Best Lensing Configurations 

In addition to identifying the general trends of lensing 
cross section with various parameters, we also examine 
how the mass configurations resulting in the very best 
(j/^i values behave. The results of Figures [8] and [9] show 
considerable scatter about the mean relation, so there 
may be parameters that can cause particular mass con- 
figurations to scatter to a high cross section. These mass 
configurations, though likely very rare in the universe, 
could potentially be the very best lines of sight for lens- 
ing high-redshift galaxies. 

We examine the frequency of 2-halo lines of sight that 
are ranked in the top 10% when sorted by af^ as a func- 
tion of various parameters. The values of afj, correspond- 
ing to the top 10% decile are given in Table [1] In Fig- 
ure [TOl we plot the fraction of trials that fall within the 



top 10% in cr^ in each bin for both projected separa- 
tion and mass ratio. This plot represents the excess (or 
deficit) of high af^ mass configurations in each parame- 
ter bin. By plotting the results in this manner, we can 
see if the parameter ranges that generate the very best 
mass configurations for lensing are similar to what would 
be expected from the mean trends for the full sample. 
In addition, this visualization removes the effects of our 
priors on the parameter being plotted such that we are 
not inherently biased toward a region of that parameter 
space that was just sampled more frequently and would 
therefore tend to have a larger raw number of configu- 
rations that fell in the top 10% due to scatter. We are 
effectively examining what values of a particular param- 
eter result in the best lensing configurations regardless 
of their actual likelihood, provided that likelihood is not 
vanishingly small. 

We see that the angular separation and mass ratio de- 
pendencies are consistent with the mean results, with the 
separation bin around ^100" and mass ratio bins near 
unity containing a larger fraction of top 10% configu- 
rations. The dependence on angular separation demon- 
strates that our choice of a 3.'5 maximum separation is 
justified, as the very best mass configurations tend to 
have separations smaller than this. 

We make similar plots for the physical properties of the 
halos in the single-halo case in the top row of Figure [TTl 
Interestingly, the frequency of top 10% mass configura- 
tions peaks at low redshifts {z ^ 0.2 — 0.3) and decreases 
with increasing redshift. This is in contrast with the 
mean trend for the entire sample, which shows lower af^ 
at low redshift with a subsequent rise and plateau at 
higher redshift. Revisiting the upper left panel of Fig- 
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Fig. 10. — Fraction of 2-halo beams in each separation (left) and mass ratio (right) bin ranked in the top 10% of the entire sample by cr^. 
The dotted line represents a frequency of 10%, which would be expected if there were no excess of top ct^ configurations in a particular 
bin. These relations match what would be expected based on the mean relations — there is an optimal nonzero angular separation where 
the best mass configurations lie, and mass ratios closer to unity generally produce good lensing configurations. The error bars are larger in 
bins with fewer samples. 



urc [9l there is a subset of trials with the very highest 
cross sections located at low redshifts, but there is also 
a large subset of trials in these bins that have very low 
af^ values. What is happening is that at these low red- 
shifts, the NFW scale radius, rg, becomes large in an- 
gular size, even though its physical scale radius may be 
comparable to higher redshift halos of similar mass. This 
implies a decrease in the surface density, which decreases 
the lensing strength of the halo. In addition, the critical 
density, Sc, is larger at lower redshift, further decreasing 
the lensing strength. However, halos that have an unusu- 
ally high concentration at these redshifts will have a large 
surface density over a large angular area, creating a favor- 
able lens. Indeed, an inspection of the high-cross section 
trials at low redshift reveals a disproportionate amount 
of halos with high concentration. In general, a halo at 
low-redshift will have an unfavorable lensing geometry 
compared to a similar halo at higher redshift, but in cer- 
tain rare cases, the concentration can be high enough 
that the halo becomes a very good lens. An example is 
Abell 1689, which has a redshift oi z = 0.18, but has an 
unusually high concentration co mpared to expectations 
from cosmological s imulations ([B roadhurst et al .l 120051: 
Halkola ct al. 200d: iZekser et aL 2006; Limous in et alj 
2007; Umetsu & Broa dhurst! " l2008t iCoe et al.n2QTor 
Sercno & Umetsu 2011), making it one of the best lens- 
ing clusters known. 

The other physical parameters that we found relevant, 
namely the projected concentration and angle between 
the major axis and the line of sight, behave in the way ex- 
pected from the mean trends — higher concentration and 
better alignment with the LOS results in better lenses, 
and thus a higher chance of attaining a very high cross 
section. 

We make similar plots in the bottom row of Figure [TTl 
but we now plot the relative fractions for the 2-halo 



beams as a function of the physical parameters of both 
the more massive and less massive halo in each trial. If 
the properties of the more massive halo dominate the 
effect on the lensing properties, as was suggested by Fig- 
ure [HI then we would expect to see similar trends with 
the properties of the more massive halo, but essentially 
no trend with those of the less massive halo. This is true 
for the concentration and for w, as the trends with the 
more massive halo are similar to those in the single-halo 
case, and the trends with the less massive halo are flat. 
However, there is a trend with redshift of the less mas- 
sive halo that is qualitatively similar to that with both 
the redshift of the high mass halo and of the single-halo 
case. This result arises at least in part from the shape 
and evolution of the halo mass function. In roughly 3/4 
of our 2-halo trials, the more massive halo in a pair is at 
lower redshift than the less massive one. In other words, 
when the less massive halo is at a low redshift, it is more 
likely for the other, more massive halo, to also be at low 
redshift. Furthermore, such configurations tend to have 
less extreme mass ratios. Since low redshift halos show 
a higher frequency of top 10% a^ values, it appears that 
the redshift of the less massive halo is having an effect, 
but in fact, this trend is primarily an artifact of the mass 
function and our requirement that the lines of sight have 
a fixed total mass rather than a real lensing effect. 

The dominance of a single massive halo in most of the 
2-halo case initially appears to be inconsistent with the 
trend showing that more equal mass ratios are more likely 
to result in a good configuration. In fact, the fractions 
shown in the right panel of Figure [10] indicate that there 
are still many excellent mass configurations that do not 
follow the trend, especially considering that the bins rep- 
resenting the most extreme mass ratios contain far more 
trials compared to the bins of more equal mass ratio. The 
point is that none of these parameters are required to be 
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Fig. 11. — Top: Fraction of single-halo beams in each redshift (left), concentration (center), and orientation angle (right) bin ranked in 
the top 10% in ai^. Bottom: Fraction of 2-halo beams in each bin of the same parameters for the more massive (black solid line) and less 
massive (red dashed line) halo in each trial that is ranked in the top 10% in cr^. The dotted line represents a frequency of 10%, which would 
be expected if there were no excess of top ct^ configurations in a particular bin. This visualization eliminates the effect of the frequency of 
certain mass configurations in the parameter being plotted on the x-axis. The concentration and cu relations for the single halo and more 
massive halo in the 2-halo case match what would be expected from the mean relations in Figure [9] Low redshift bins (z ~ 0.2) contain a 
larger fraction of top 10% mass configurations, despite the mean trend tending to disfavor low redshift. The physical properties of the less 
massive halo in the 2-halo case are generally unimportant since there is usually a single dominant halo, even in the top 10% trials. There 
is a preference for the lower mass halo to be at low redshift, but this arises at least in part from a selection effect where such configurations 
are also likely to have the higher mass halo at low redshift, and the solid line indicates that having the higher mass halo at low redshift is 
more likely to result in a top 10% Cfi configuration. 



optimized in order to produce a good mass configuration. 
Given a fixed total mass, there are certainly parameter 
values that increase the likelihood of this occurring, but 
there are multiple ways and combinations of parameters 
that can result in a large cross section. 

5. DISCUSSION 

5.1. Lensing Interactions Between Halos 

Our results show that there is a gain in the lensing 
cross-section when a massive halo is split into two or 
more cluster-scale halos projected close on the sky. How- 
ever, the results presented in ij 14. II indicate that it is not 
simply the lensing properties of the individual halos that 
sum up to produce a more powerful lens than a single 
halo of equivalent total mass. In fact, the values shown 
in Table [1] and Figure [7] demonstrate the opposite ef- 
fect. Instead, the gain in cross-section for the multi-halo 
beams over the single-halo case of equivalent mass arises 
from interactions between the lens potentials of the ha- 
los, even if the halos are not physically interacting. The 
results presented in this paper are specific to our chosen 
fiducial mass and range of halo separations, but some 
general statements can be inferred. 

There is an optimal separation (~ 100" for our cho- 
sen fiducial total mass) between halos that maximizes 



the cross-section (Figure [TU)) . One can understand this 
conceptually as follows. Much of the contribution to cr^ 
comes from regions of the source plane that are outside 
the strong lensing caustics but still significantly magni- 
fied (see Figures [5] and HI). Outside the caustics, the mag- 
nification produced by a given halo generally decreases 
with distance from the halo. When two or more halos 
are placed in proximity to each other on the sky (even 
if they are at different redshifts) , their magnification re- 
gions in the source plane can overlap and create areas 
where no halo individually exceeds our chosen magni- 
fication threshold but the joint lensing effects boost /i 
above the threshold. While this notion of overlap is use- 
ful heuristically, recall from § 12.11 that the net magnifi- 
cation is not simply the sum of the magnifications from 
individual halos. As indicated by Equations HHTUl we 
need to consider not only the overlap of the halos' sur- 
face density profiles (the k terms), but also the combined 
effects from shear (the 7 terms). In other words, the 
"lensing interaction" between the halos is mediated by 
both convergence and shear. 

As the halos' angular separation decreases, the boost 
from the lensing interaction begins to increase. At some 
point, however, the separation becomes small enough 
that the caustics begin to merge on the sky, which may 
increase the area where the magnification is very large 
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(e.g., fj, > 10-100) at the cost of reducing regions where 
the magnification is more modest but still significant 
(e.g., n ^ 3-10). In other words, finite separations 
can maximize the area that produces modest magnifi- 
cation, whereas smaller separations could be more rele- 
vant if we were interested in larger magnifications (see 
Appendix r^ . 

A related effect explains why having multiple halos 
helps. As a beam containing a fixed total mass is split 
into multiple halos, each halo's average mass decreases 
and so its magnification region shrinks. At the same 
time, the average projected separation between halos also 
decreases (when the beam size is fixed). For the beam 
sizes and masses we have considered, the decrease in a^^ 
from the dilution of mass is more than offset by the in- 
crease in af^ from the proximity of halos, so the cross 
sections tend to increase when we consider beams with 
more halos but the same total mass. The etendue of the 
mass configuration, defined as magnification over a large 
area, is increased as the mass is dispersed, provided that 
it is not spread out so thin that regions in between the 
halos drop below the magnification of interest. 

We have examined multi-halo trials where we have 
changed the total integrated mass, and while the optimal 
offset may change, our general conclusion that splitting 
the mass into multiple halos holds true, provided that 
the typical angular separation between halos is roughly 
comparable to the size of their region of infiuence. 

5.2. Additional Sources of Uncertainty 

We have generated mass configurations that are phys- 
ically motivated by the relative frequencies and proper- 
ties of cluster-scale halos as determined from simulations, 
with the exception that we have constrained these lines 
of sight to have a large total mass to represent power- 
ful lensing configurations. Although our treatment of 
these massive halos accounts for much of the diversity 
seen among such halos with regard to properties such 
as redshift and mass distribution, concentration, and 
shape/orientation, among other things, there are some 
details that we have not considered. 

While we allow for multiple dark matter halos at sim- 
ilar redshift in our models, we do not account for sub- 
structure in our halos. This includes the effects of galax- 
ies and higher-order details in the halo models (e.g. 
misaligned isodensity co ntou rs, mass-shape correlations; 
Hchneider et al. 20111). IMbert et al. (2008) find that 
the lensing cross sections for halos more massive than 
10^"* M0 increases by at most 40% for a source redshift of 
z — 5.7 and is more significant for lower masses when the 
effects of galaxies are included, although their definition 
of cross section is somewhat different from ours. The 
galaxies will affect both the single-halo and multi-halo 
cases, and in fact could have a larger impact on the multi- 
halo fields where the typical halo mass is less than the 
fiducial mass, so our assumption of no substructure due 
to a galaxy component may be conservative. The effects 
of baryonic cooling can increase the lensing cross section 
of these halos (R,ozo et al. 2008 ), but this should affect 
both single-halo and multi-halo lines of sight in a similar 
fashion. We also assume that our halos follow an NFW 
density profile, and deviations from NFW could alter the 
calculated cross sections. However, as noted earlier, this 
assumption makes our analysis compatible with previ- 



ous work and is suggest ed by some observational results 
(Kne ib et al.|[20 03: Kat gert et an[200l [Comerford et al.l 
l2006tlShuetal.ii2008.) . 

Although we allow for the possibility of halos being 
at the same redshift, we do not account for physical in- 
teractions between halos that will alter the dark mat- 
ter distribution, such as mutual tidal effects, density en- 
hanceme nts in the baryoni c components (e .g., the Bullet 
Clust er (jClowe et al.ll2006D or Abeh 2744 (jMerten et al.l 
I201in ). and multiple subclusters sharing a common dark 
matter halo. These effects may lead to an overall un- 
derestimate of the lensing cross section, but should not 
preferentially affect a given multi-halo configuration. 

We also do not account for smaller uncorrelated struc- 
tures along the line of sight. These low-mass structures, 
which include field galaxies and smaller groups, should 
have small effects on the lensing cross section of the main 
cluster (s) in these fields if the results of Figure |S] are ex- 
trapolated to lower mass ratios. We do not expect such 
structures to be more prevalent in beams with fewer ha- 
los, so it will not affect our general conclusions about 
multi-halo fields. We have required that the halos in our 
analysis have a minimum mass of 10^'^Mq, but less mas- 
sive halos down to the group scale (~ lO^^M©) are more 
common and are more likely to be found in multi-halo 
configurations. We run a test where we set the minimum 
halo mass in the 2-halo case to 5 x IO^^Mq and find that 
(Tf^ is still larger for these cases on average than for the 
single-halo case. This also can be seen from the trend 
in mass ratio, where a larger difference between the two 
halo masses results in a smaller a^. 

5.3. Implications for Detection of High-Redshift 
Sources 

Gravitational lensing is effectively a mapping from a 
given area in the source plane to an area in the image 
plane, which aids in the detection of faint sources by in- 
creasing the solid angle subtended with no loss in surface 
brightness. However, this magnification comes at the 
cost of a decreased area being probed in the source plane. 
As a result, the gains achieved from lensing in terms of 
detecting high-redshift galaxies depends very strongly on 
the shape of the source galaxy luminosity function (LF) . 
The LF i s typically assumed to be of the form first pro- 
posed bv lSchechteii (|1976D . which is parameterized by a 
characteristic magnitude M*, a faint-end slope a, and a 
normalization factor (jj* representing the number density 
at the characteristic luminosity. For galaxies fainter than 
M*, the LF is approximately a power law with slope a. 
For galaxies brighter than M* , the LF declines exponen- 
tially with increasing luminosity. 

Since the loss in area probed is inversely proportional 
to the lensing magnification, lensing increases the like- 
lihood of detecting a galaxy at a given source redshift 
if the LF at that galaxy's luminosity is declining more 
steeply than a power law of index —2 (not accounting 
for finite source size effects; see Amnions et al. in prepa- 
ration). Therefore, lensing is mainly bene ficial for detect- 
ing the most lum i nous sources ( Broadhu rst et al.l 119951 : 
ISantos et al.ll2004lBouwens et al. 2009), unless the faint- 
end slope is steep. The shape of the LF is not well con- 
strained at z > 7, so it is not immediately obvious how 
(7^ would translate to an enhancement in the detectabil- 
ity of high-redshift sources. 
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We have presented our results in terms of lensing cross- 
section or etendue. Observational studies of the high- 
redshift universe have a variety of science objectives that 
may require different observing strategies in terms of 
depth versus area or volume. For example, given a fixed 
amount of observing time, constraining the LF faint-end 
slope requires a very deep observation, while constraints 
on the bright-end are optimized with shallower obser- 
vations over a wider area. Lensing will detect fainter 
sources than a similar blank field observation within the 
same integration time. Here, we examine the specific 
case of detecting 2; ~ 10 galaxies and the extent to which 
multi-halo beams provide an increase in the frequency 
and depth of detected sources relative to blank fields and 
single lensing clusters. 

The detailed calculation depends on several variables, 
including the source size and ellipticity distribution, the 
instrument PSF, the depth and field of view of the ob- 
serving program, and the source luminosity function. 
The only observational cons traint on the LF at z ~ 10 
is based on ~ 0.8 detections (jBouwens et al.|[2Qllal) with 
ultra-deep HUDF observations using ^ffST/WFCS. Be- 
cause of these weak constraints, a calculation of the num- 
ber of z ~ 10 detections that would be expected in fields 
containing multiple cluster-scale halos is extremely un- 
certain. However, we can make some simplifying assump- 
tions and estimate the number of detections at different 
luminosities relative to a blank field of equivalent size. 

The enhanced detectability of high-redshift sources 
due to lensing does not scale linearly with magnifi- 
cation due to the effects of finite source size (Am- 
mons et al. in preparation). z ^ 1 sources 
are marg inally resolved in high- r esolution space - based 
imaging ([B arkana fc Loeb' 2000; 'Fergus on et al.l 120041 : 
IBouwense t al. 2008; Ocsch et al. 2010a|) and become 
smaller at higher redshifts if a scaling of i?e// oc (1 + 
2)"^-^ is assumed (|Bouwens et al.|[2006l ). For unresolved 
sources, the effects of lensing magnification can increase 
their solid angles to sizes larger than the detector PSF. 
For a lensed source resolved in one, two, or neither direc- 
tions on the sky, the fractional increase in its S/N com- 
pared to an unlensed source, or "signal-to-noise boost" 
(SNB), can scale differently with magnification in differ- 
ent regimes. This effect must be considered when calcu- 
lating the detectability of high-redshift sources. 

We focus on the case of a hypothetical space-based 
search for high-redshift galaxies and assume a PSF full- 
width at half-maximum of 0.15", approximately the size 
of the i/ST/WFCS PSF in the Rxmw filter. We as- 
sume a population of sources with a Gaussian shape 
at 2: = 10 with a normal distribution of effective radii 
with a mean of 0.12" and standard deviation of 0.04". 
This angular size distribution is determined by app lying 
a size evolution of (1 -I- z)~^-^ (jBouwens et al.|[200Q) to 
the half-light radius distribution of z ~ 7 sources, which 
is 0.75 ± 0.23 kpc for gal axies with luminosities in the 
range 0.3L* < L < L* (jOesch et a l. 2010a). We as- 
sume a normal distribution of intrinsic major-to-minor 
axis ratios with a mean of 1.2 and a standard deviation 
of 0.3, which can be folded into th e SNB calculation using 
the formalism in iKeetonI ()2001b( ) . The effect of intrinsic 
shape on the detectability of a source is much smaller 
than the effect of source size. 



comparable de pth to that of lBouwens et al.) (29.8 mag in 
WFC3 jJifinM-": l2011aD relative to a similar observation in 



The field of view and observational depth will also af- 
fect the number of detected sources. We calculate the 
number of sources in our source plane area that have at 
least one image detected in a 5' x 5' area mos aic survey of 
e pth to tl 

sinni) 

a blank field of equivalent area. We assume the z - ^ 10 lu - 
minosity function constraints of lBouwens et all ()2011al) , 
who chose fixed values of (jf = 0.0012 Mpc~'^ and 
a — —1.74 and derived a constraint of M* = — 18.3± 0.5 
mag. 

In Figure [T2I we plot the number of detections per unit 
magnitude as a function of source luminosity for both 
blank fields and our model mass configurations. The 
blank field observation can detect more objects at the 
bright end due to the larger volume probed in the source 
plane. The faint end, on the other hand, is beyond its 
reach, but can be detected in the lensing fields. The rel- 
ative tradeoff between these two detected source popula- 
tions favors faint end detections as the number of halos 
increases, reflecting the gain in the lensing cross section 
for these mass configurations. Therefore, the mass con- 
figurations with more halos, in addition to having larger 
cross sections, also probe larger source plane volumes at 
the faint end of the LF. The effects of changing various 
halo properties are generally smaller than the effect of 
changing the number of halos. 

As a result of the large uncertainty in M* (±0.5 mag), 
the uncertainties in the curves in Figure [12] due to error 
in the LF are comparable to the differences among the 
various curves. For example, in the particular case of 
M* — —18.3, the 1-halo and 5-halo beams result in a 
factor of ^^3.5 and ~4 increase, respectively, in the num- 
ber of detections fainter than M* over the blank field. 
However, changing M* by la to —17.8 results in a factor 
of ~23 and ^^37 increase, respectively, over the number 
of faint blank field detections. These large deviations re- 
sult from the fact that there are very few faint-end detec- 
tions expected from blank field observations. Better con- 
straints on both the luminosity and size distributions of 
high-redshift sources are clearly needed, and can in fact 
be provided by observations of such multi-halo beams. 

We caution against using these numbers as more than 
rough estimates due to the large LF uncertainties and 
the strong dependence on the details of the observing 
program. For example, ground-based surveys, while still 
inferior to high-resolution space-based imaging, will ben- 
efit more from lensing due to the larger PSFs of current 
ground-based detectors. Large uncertainties arise be- 
cause the assumed depth of the observation only reaches 
a luminosity that is still in the exponential part of the 
LF, so a change in M* leads to a large change in the 
gain due to lensing. Nevertheless, the improvement due 
to splitting the beam's mass into more than one halo is 
clear. 

5.4. Likelihood of Multiple Clusters in Projection 

Although we have demonstrated that lines of sight con- 
taining multiple cluster-sized halos nearby in projection 
can be better at lensing high-redshift galaxies into de- 
tectability than those containing a single cluster, it is rea- 
sonable to consider the likelihood of such fields existing in 
the real universe. A handful of serendipitous discoveries 
of cluster-cluster lenses have been reported in the liter- 
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Fig. 12. — Number of 2 = 10 detections per unit magnitude for hypothetical 5' X 5' space-based blank field (black lines) and lensing field 
(colored lines) observations. The luminosity function of Bouwcns ct al. (2011a) is assumed. The blank field observation can detect more 
objects at the bright end due to the larger volume probed in the source plane, while the faint end is only detectable in the lensing fields. 
The tradeoff between the bright and faint end tends more toward faint end detections as the number of halos increases, refiecting the gain 
in the cross section. The mass configurations with more halos, in addition to having larger cross sections, also probe larger source plane 
volumes at magnifications needed to detect the faint end of the LF, resulting in more detections of objects fainter than M*. 



ature (jBlakesled [200lL lAthreva et all I200I iZitrin et all 
I2012D . albeit with smaller total integrated masses than 
our simulated fields, so their existence is not simply a 
theoretical construct. Here, we perform a simple cal- 
culation using our assumed halo mass function to esti- 
mate the number of fields in the universe that contain 
two cluster-scale halos with a total integrated mass of at 
least 2 X IO^^Mq. 

We begin by performing the calculation of 2-halo 
beams in general terms, then examine specific cases that 
are relevant in the context of our simulations. We focus 
on the case of two halos for simplicity, requiring that each 
halo has at least a minimum mass, Mmin, and lies within 
the redshift range Zmin "£ zi < Z2 < z„iax- Multi-cluster 
lines of sight will be strongly biased toward low masses 
due to the declining mass function, so we also require 
that the two masses sum to some minimum total mass, 
Mtot- We define the quantity 



n(M , z) 



dN 
dMdV 



= in 



H{z) dN 
c dMdz 



(12) 



to be the number of halos with mass M at redshift z 
across the entire sky, where H{z) is the Hubble parame- 
ter at redshift z. Using the halo mass function to deter- 
mine dN/dMdz, we then calculate the number density 
of foreground halos, n(Mi,zi), and the number density 
of background halos, n{M2,Z2), by evaluating the mass 
function at the appropriate masses and redshifts. We 
account for the maximum separation range of the two 
halos by requiring that the second halo lies within a cir- 
cle of angular radius Rsep- This is a fixed quantity for 
all redshifts, and therefore is a constant that factors out 
of the integral. If we had chosen a geometry-dependent 
maximum separation (e.g., the Einstein radius of the 
foreground lens), we would need to include the solid- 
angle term when integrating over Z2 ■ Our restrictions on 



both the minimum halo mass and total mass restricts the 
minimum value of M2 to either Mmin or (Mtot — Mi), 
whichever is greater. With these assumptions, the total 
number of 2-halo lines of sight that have a mass of at 
least Mtot is given by 



N = 



Att ster 



dMi 



'M„,i„ 



/ dM2 

'inax(A/„ji„, A/tot -i\/i) J zi 



dzi 

n{Mi, zi)n{M2, Z2)dz2- 
(13) 



With this general calculation, we can then examine 
specific cases. We look at the 2-halo case analyzed in our 
Monte Carlo trials by setting M„„ = lO^Mg, Mtot = 

2 X IO^^Mq, Zmin = 0.1, Zmax = 0.9, and Rsep = 3.'5. 

Based on this calculation, we expect there to be ^10 
such lines of sight in the sky. This indicates that there 
exist ^1 field similar to the top 10% of just the 2-halo 
beams analyzed in S I4.31 so it is reasonable to expect that 
such powerful lensing beams do exist. The frequency of 
multi-halo lines of sight also increases dramatically if our 
requirements on total mass or minimum halo mass are 
relaxed. 

The calculation is very sensitive to our assumed cos- 
mology, particularly as- The uncertainty in erg from 
the seven-year WMAP results suggests a factor of ^ 2 
uncertainty. This simplistic calculation is likely conser- 
vative, as it does not account for mass due to smaller 
structures along the line of sight, the possibility of more 
than two cluster-scale halos in projection, or the effects 
of large-scale structure that could lead to correlations 
among mass peaks. Accounting for these effects could in- 
crease the predicted total number of massive multi-halo 
fields. A more detailed analysis of the frequency of these 
high-mass lines of sight using the Millenium Simulation 
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will be presented in French et al. (in preparation). Ad- 
ditionally, we have identified several high-mass lines of 
sight containing multiple structures from spectroscopic 
observations of selected SDSS fields. An ongoing spec- 
troscopic and photometric observing program has discov- 
ered and confirmed at least two lines of sight that contain 
multiple halos each with dynamical masses greater than 
~ 10^'* Mq as inferred from the relative velocities of the 
cluster galaxies. The total integrated masses in these 
fields are estimated to be > 3 x lO^^M©. The details 
of these observations will be presented in Amnions et al. 
(in preparation). 

6. CONCLUSIONS 

We have compared the lensing properties of lines of 
sight containing chance alignments of multiple cluster- 
scale halos on the sky to single-halo lenses, isolating the 
effects of more projected halos and the detailed prop- 
erties of those halos. We examine the implications for 
magnifying very high-redshift (z ^ 10) galaxies. Our 
main conclusions are as follows. 

• Lines of sight containing multiple cluster-scale ha- 
los have larger cross-sections for lensing high- 
redshift galaxies above a given magnification 
threshold than single halos of equivalent total in- 
tegrated mass. This effect is driven by interactions 
among the lens potentials of the halos, arising from 
non-linear combinations of both their convergence 
and shear profiles. The gain in the lensing strength 
due to the presence of multiple halos is maximized 
at some finite projected separation. For a beam 
with a total mass of 2 x lO^^M©, this separation is 
~100". At larger separations, the configuration ap- 
proaches the limit where the halos are independent 
fields. At smaller separations, the mass distribu- 
tion is inefficient at maximizing the field's etendue 
(the area containing large magnifications), result- 
ing in a lower cross section. 

• For a line of sight containing the same total inte- 
grated mass split between two cluster-scale halos, 
the gain in lensing cross section is larger in those 
rare cases where the halo masses are comparable. 

• The physical properties of a halo that drive an in- 
crease in the lensing cross section are a high con- 
centration, the major axis aligned with the line 
of sight, and an intermediate redshift (> 0.4), 
which optimizes the scaled surface density of the 
halo. The triaxiality and projected ellipticity of 
the halo, and the relative orientations of multiple 
halos, are relatively unimportant. For beams with 
two cluster-scale halos, one halo typically domi- 
nates the mass due to the shape of the halo mass 
function. In such cases, the properties of the dom- 
inant halo determine the lensing characteristics of 
the field. One implication of this result is that 
analyses of known single-halo lensing beams, which 
rarely account for other structures along the line of 
sight, are unlikely to have a second massive halo in 
the beam that would significantly alter the overall 



lensing cross section, but errors in the derived 2-D 
magnification map may still be underestimated. 

• The same factors that increase the lensing cross 
section in typical 2-halo beams generally produce 
the best (top 10%) of 2-halo beams. One excep- 
tion is that halos at low redshifts that scatter to 
large concentrations, resulting in an increased sur- 
face density over a large area, can have unusually 
large lensing cross sections. This is opposite to 
the overall trend for halos at low redshift to have 
smaller cross sections. 

• We have isolated the effects of mass configuration 
— multiple halos and the properties of those ha- 
los — while controlling for total integrated mass in 
the beam. These effects are non-negligible and can 
be comparable to increasing the total mass in the 
beam from - lO^^M© to ~ 3 x IO^^Mq. Actual 
lines of sight in the universe containing multiple ha- 
los have the additional benefit that they are likely 
to contain more total mass than those with just a 
single massive cluster (French et al. in prepara- 
tion). 

• We have conservatively estimated the number of 
lines of sight that contain two cluster-scale halos 
within 0.1 < z < 0.9 as ~10 in the sky. A more 
detailed constraint using cosmological simulations 
will be presented in a future work (French et al. in 
preparation). We have identified several fields in 
the SDSS with large integrated masses and multiple 
cluster-scale halos in projection. Detailed observa- 
tions of these fields will be presented in Ammons 
et al. (in preparation). 

The results of our analysis can be used to identify ideal 
lines of sight in the real universe for magnifying very 
high-redshift galaxies into detectability once the masses, 
redshifts, centroids, and certain physical properties (e.g., 
concentration) of multiple cluster-scale halos are con- 
strained from observational data. This work shows that 
focusing on reducing the uncertainties in these parame- 
ters will most limit the errors in the resulting magnifica- 
tion maps, which in turn will improve constraints on the 
luminosity function of lensed sources at z ~ 10 identified 
in these fields. 
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Fig. 13. — Left: ct^ relative to the 1-halo case as a function of the threshold magnification, fit at which a^, is defined for each of the 
beam configurations (top). We make a similar plot, including only the top 10% of beams ranked by aij,^fj,^ (bottom). The gain in cross 
section for the multi-halo fields over the single-halo fields is always pos itiv e on average, regardless of our choice of magnification threshold 
for defining tr^t in this range. Right: Same as the left panel of Figure [Tol but for cr^ defined to include magnifications above 3 (black), 5 
(red), and 10 (blue). The peak shifts to smaller separations as fit increases, reflecting the fact that larger projected surface densities are 
needed to produce regions of at least the threshold magnification value. 



A. 



APPENDIX 
CHOICE OF (jf, MAGNIFICATION THRESHOLD 



Our chosen threshold of /it = 3 for the brightest image of a given source position is somewhat arbitrary, ahhough 
there is evidence that the number of detections in lensing fields as a function of magnification broadly peaks around 
yU ~ 3 (Ammons et al. in preparation). However, one could imagine that for the purposes of magnifying high-redshift 
galaxies, we may want a different magnification threshold. At z ~ 10, realistic observational programs with HST or 
a future facility will have little or no ability to probe significantly fainter than M* . In the exponential part of the 
luminosity function, the gain in source density outweighs the loss of volume probed for increasing lensing magnification, 
so higher magnification may be desirable. 

We reevaluate a^ for different magnification thresholds, ranging from /zt = 3 to /it = 20, and plot the mean values 
relative to the single- halo case in the top panel of the left plot of Figure [T31 We can see that the choice of magnification 
threshold within this range does not change our general result that the multi-halo fields have larger cross-sections than 
the single-halo case, although the fractional gain over the single-halo case does vary. In the bottom panel, we plot 
the same quantity for just the top 10% fields. We do not examine higher magnification thresholds because the chosen 
resolution of our source plane grid starts to become problematic, and it would require a large amount of computation 
time to reevaluate these configurations on a more finely sampled grid. Likewise, smaller magnification thresholds 
(1 < /^t ^ 3) require a calculation over a larger source plane area and are starting to become less scientifically 
interesting since they are approaching the case of no lensing magnification. The results appear unlikely to change at 
low redshift anyway given the nature of the trends. 

Another possible consequence of changing nt is that the the effects of various parameters on a^, such as the ones 
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discussed in § 14.21 may change. In general, increasing this threshold will reduce the normalization of cr^, but it is 
uncertain whether its behavior as a function of these parameters will change. Performing the simple mass scaling test 
shown in Figure [T] for different choice of fj,t only appears to change the overall normalization, but the nature of the 
trend remains the same. What about the other physical parameters that were explicitly tested in § 14.21 ? 

We test the effects of changing /if for various physical parameters and find that they are qualitatively similar with 
the exception of angular separation. In the right panel of Figure [TH we replot the left panel of Figure [101 but with the 
results for /it = 5 and /it — 10 as well. We plot this figure rather than one analogous to the upper left panel of Figure[S] 
to remove the effect of overall normalization, and because the trends with separation that applies to the top 10% of 
beams are qualitatively similar to the mean results. As we increase /Zt, the optimal angular separation decreases. This 
makes sense because in order for interactions between the lens potentials to produce extended magnifications regions 
above an increasing threshold, the interactions between the convergence and shear terms of the different halos need to 
be stronger. This indicates that the halos need to be projected closer together since their projected density and shear 
profiles drop off with increasing radius. 

B. PROJECTION EFFECTS ON HALO PROPERTIES 

The lensing properties of a given mass distribution are determined by its projected mass density in the plane of the 
sky. In our analysis, we assume NFW halos with a triaxial mass distribution, which we need to translate into projected 
qua ntities. Here, we describe ou r method for char acterizing the projecte d mass distribution of a triaxial ellipsoid ( see 
also lde Zeeuw fc Fr anx 1989; O guri et al . 2003; va n de Ven et al...2009:.Morandi et al.l[201lL iFteroz fc HobsonllMa ). 

For a spherical NFW halo, the density profile is given by 

[r/rs)(l + r/rsY 

where r^. is the scale radius where the slope of the mass profile matches that of an isothermal sphere, and ps is the 
characteristic density. When projected onto the plane of the sky, the surface density profile is given by 



^q{r)= [ p Ur'^ + zA dz 

/•OO 

= 2 / p{R cosh u) R cosh udu 
Jo 



(B2) 



with the change of variables z = i?sinhu. Since p{r) depends on its argument only through r/rg, we know that T,q(R) 
depends on R only through the combination R/r^. 

For a triaxial NFW halo, we specify its major, intermediate, and minor semi-axes by a, b, and c, respectively. To 
simplify the notation in this section, we rescale these semi-axis ratios by r^ such that a = a/vs , b = b/r^ , and c = c/rs 
are dimensionless quantities. Let X = {X, Y, Z) be coordinates in the halo's principal axis frame, and x = {x, y, z) 
be coordinates in the observer's frame such that the z-axis is along the line of sight. The two coordinate systems are 
related by a rotation matrix, 

X = Ax. (B3) 

We then define 

2 X2 y2 ^2 

"^ = — + -^ + — ' (B4) 

a c 

where m has dimensions of length and is analogous to r in the sph erical case. We can define a triaxial NFW halo by 
substituting m into Equation IBll to get p{m). Using Equation IB31 we can write m^ as a quadratic in {x,y,z): 

m^^fz^+gz + h, (B5) 

where /, g, and h are functions of x, y, and the angles used to define the halo orientation. Explicit expressions for 
these quantities are given in § IB. II and IB. 21 We rewrite Equation IB5I as 

m2=e'+C', (B6) 

where 

C-/^/^(- + ^), (B7) 



and 




(B8) 



Now when we do the projection, we can change integration variables from z to ^. All of the {x, y) dependence is in ^, 
and C,'^ is quadratic in {x,y) such that there is some rotated coordinate system {x',y') where 
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This means that any function of C wiU have elhptical symmetry. The quantities a! and h' are the semi-major and 
semi-minor axes (again scaled by r^ so they are dimensionless quantities), respectively, of the projected ellipse. The 
frames are rotated by some angle f/'j so 



V 



cos '0 sin -0 
— sin-;/' COS';/' 



We can rearrange Equation IB8I and write it as 



S — ^xx^ ~r ^xy-^y ~r ^yyU 



(BIO) 



(Bll) 



where d^x, d^y, 



^ and dyy are the coefficients of the x^, x y, an d y^ terms, respectively. Equating Equation IBllI with 

Equation IB9I and substituting for x' and y' via Equation lBlOl we can solve for 0, a', and h': 



lb = — arctan 
2 



^xy 



where 



6' = 

F = 
G = 



F-G' 



F + G' 



J-yyi 



y [d-xx dyy) + d^y. 



(B12) 
(B13) 
(B14) 

(B15) 
(B16) 



We can implement one of two equivalent approaches using this basic framework. We can specify the orientation of 
the halo in terms of two viewing angles and apply a random rotation on the sky at the end. Alternatively, we can 
specify the halo orientation in terms of three Euler angles which will naturally produce a random projected orientation. 

B.l. Viewing Angle Approach 

Three angles are needed to fully specify the orientation of a halo in three dimensions. In the viewing angle approach, 
we only specify two angles, so the projected mass distribution must be rotated by a random angle at the end to achieve 
different position angles. In this approach, we define 6 and 4> to be the polar and azimut hal angles, respectively, of 
the ellipsoid relative to the observer's frame. Following the formalism in lOguri et al.l ()2003l ). the rotation matrix A is 
given by 



A = 



The coefficients in Equation IB5I are given by 



— sin0 


— cos cos (j) 


sin 9 cos 


coscf) 


— cos 9 sin (j) 


sin 9 sin 





sin 9 


cos6' 



/ = sin^ 



cos 



sm 



62 



cos 



sin sin 20 [ —r 

,„2 



h 



52 



sm 29{ —- 



€05^9 



cos 



sm 



62 



y, 



a? 62 J c 

Proceeding with these definitions, the coefiicients in Eg nation IB 111 are given by 

1 



y + cos 9 sin 20 



1 
52 



xy. 



^vv 



D 

hi' 

1 



[a^ cos^ 9 cos^ ( 



6^ cos^ 9 sin^ 1 



c^ sin^ ( 



2 • 2 
a sm ( 



b^ cos^ ( 



dxy = — {b ~ a ) cos sin 20, 



D = a'^b'^ cos^ 



2 • 2 
c sm I 



(a^ sin^ ( 



i -I- 6^ cos^ 



(B17) 

(B18) 
(B19) 
(B20) 

(B21) 

(B22) 

(B23) 
(B24) 



The angle "0 is not relevant in this approach because we have only specified two angles, which is not sufficiently general 
to fully characterize the halo orientation. Once the axis ratios of the projected ellipse have been determined, the ellipse 
can be rotated by a random angle to achieve different position angles. 
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B.2. Euler Angle Approach 

Alternatively, we can specify the orientation of the halo by three Euler angles, a, /?, and 7. These angles are related 
to the viewing angles 9 and cj) by the transformation 



a — 0, 

/3 = -i 



2' 



(B25) 
(B26) 
(B27) 



In this formulation, the rotation matrix becomes 



cos a cos 7 — sin a cos /3 sin 7 cos a cos /? sin 7 + sin a cos 7 sin /? sin 7 

— sin a cos /3 cos 7 — cos a sin 7 cos a cos /3 cos 7 — sin a sin 7 sin /3 cos 7 
sin a sin /3 — cos a sin /3 cos/3 

The quantities /, g, and h can then be expressed in terms of the three Euler angles as 

cos^^ 



f SIT? 7 cos^ 7 



9-- 



+ 



1 l\ ■/,•„ . ■r,/,/^sm7 cos 7 
^r — TTT cos a sm p sm 27 — sm a sm 2p — ^ 1 -7: — 



1 



1 



1 



-^ 7: \ sin a sin /3 sin 27 + cos a sin 2/3 



sin^ 7 cos^ 7 



(cosacos7 — sinQ;cos/3sin7)2 (sin a cos /3 cos 7 + cos a sin 7)^ sin2Q;sin2/3 



(cos a cos /3 sin 7 + sin a cos 7)2 (cos a cos /3 cos 7 — sin a sin 7) ^ cos2Q:sin^/3 



a2 62 

2 (cos a cos /3 sin 7 + sin a cos 7) (cos a cos 7 — sin a cos /3 sin 7) 



2(sin a cos /3 cos 7 + cos a sin 7) (sin a sin 7 — cos a cos /3 cos 7) sin 2a sin /3 

62 ^^2 



xy. 



(B28) 



(B29) 
(B30) 

(B31) 



The coefficients dxx, dxy, and dyy in this approach are 



dxx = — [(sin a cos 7 + cos a cos /3 sin 7) 2a^ + (sin a sin 7 — cos a cos /3 cos 7)^6^ + (cosasin/3)2c2] 
dyy = — [(cos a cos 7 — sinacos/3sin7)2a^ + (cos a sin 7 + sin a cos /3 cos 7)^6^ + (sinasin/3)2c2] 



1 

D 

,21,2 



sin 2a sin /3 c 



{a^ —b"^) ( - (3 + cos 2/3) sin 2a cos 27 + cos 2a cos /3 sin 27 



D = a^h^ cos2 /3 + c2 sin^ I3{a^ cos^ 7 + 6^ sin^ 7). 



(B32) 
(B33) 

(B34) 
(B35) 



The angle 1/' (Equation |BT2| is now the position angle of the projected surface mass distribution, whereas a random 
position angle was drawn independently in the viewing angle approach. Both approaches give equivalent distributions 
of projected ellipticities and position angles. 



B.3. Effects on the Density Profile 

The relevant physical quantities for determining the lensing properties of an NFW halo are the ellipticity e, the scale 
radius rj, and the scaled surface density Kg — VsPs/'^ci where Sc is the critical surface density for lensing. Assuming 
a triaxial rather than spherical halo does not change its virial radius or mass, but does change e and the projected r^ 
and Ks- The projected ellipticity is simply 



6' 
e = 1 :■ 



(B36) 
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To obtain the surface density profile of a triaxial halo, we substitute C into Equation IB2I 

S(C)=/ p{m)dz 

= /"^ X 2 / p(C cosh u) C cosh wdu 
Jo 

- /-1/2 X Eo(C), 

with the change of variables ^ = Csinhw. In other words, we can obtain the function for the surface mass density of 
the elliptical case just by plugging the elliptical radius C, into the function Eq and rescaling by f~^^^- We know that 
the position coordinates enter only through the combination C/^s- In the coordinates aligned with the major axis of 
the ellipse, this has the form 



c 



X 

a'r. 



y 

b'r. 



,-,1/2 



(B38) 



We can therefore think of the projected mass distribution having some scale radius r^, and if we choose this so a circle 
of radius r^ has the same area as the ellipse above, then we have 



Tj, — TsVa'b'. 
Since the virial radius of the halo remains the same, the projected concentration becomes 

r'. ~ sH^' 



(B39) 



(B40) 



The scaled surface density k^ is rescaled by a factor of / '^/^ as per Equation IB371 and can also be expressed as a 
combination of 3D and projected axis lengths (jvan de Ven et al.ll2009( ). 



/■ 



-1/2, 
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(B41) 
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